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The contribution of B meson decays to non-photonic electrons, which are mainly produced by 
the semi-leptonic decays of heavy flavor mesons, in p + p collisions at y^ — 200 GeV has been 
measured using azimuthal correlations between non-photonic electrons and hadrons. The extracted 
B decay contribution is approximately 50% at a transverse momentum of pt > 5 GeV/c. These 
measurements constrain the nuclear modification factor for electrons from B and D meson decays. 
The result indicates that B meson production in heavy ion collisions is also suppressed at high pr- 

PACS numbers: 25.75.-q 



The suppression of non-photonic electron yields from 
semi-leptonic decays of D and B mesons for pT up to 9 
GeV/c in central Au+Au collisions at Relativistic Heavy 
Ion Collider (RHIC) has been observed to be large [1[, 
and similar to that of light quark hadrons [2|. Due to 
the dead cone effect, heavy quarks were expected to lose 
less energy than light quarks [3[ if the dominant energy 
loss mechanism is gluon radiation [J]. Various models 
have been proposed to explain the large suppression of 
non-photonic electron yields [5|-l7[. Theoretical calcula- 
tions of the non-photonic electron suppression crucially 
depend on the B/D ratios because the amount of radia- 
tive energy loss depends on the quark mass. Measuring 
the bottom quark contribution to non-photonic electron 
yields uYp-\-p collisions is therefore important in order to 
understand the production of heavy quarks and to pro- 
vide a baseline for the energy loss measurement of heavy 
quarks in the hot and dense medium produced in central 
Au-I-Au collisions. 

In this paper, we report a determination of the rel- 
ative contribution from B decays to non-photonic elec- 
tron yields (r^ ) by measuring the azimuthal correlations 
between non-photonic electrons and charged hadrons 
[e-aon-i-h) , and between non-photonic electrons and D^ 
mesons {enon^-D^) in p + p colhsions at y/s = 200 GeV 
by the STAR experiment at RHIC. We fit the experimen- 
tal £110117-^ correlations using a combination of PYTHIA 
calculations [8[ for D and B meson decays and extract 
rs as a function of pt (2.5 < px < 9.5 GeV/c). An 
independent measurement of the rg is obtained from 
Gnon-y-D'^ Correlations, by selecting the charge combina- 
tions e~-I?"(— > K~) and e^-D (— >■ K^), which provide 
relatively pure samples of B decays and charm pairs on 
near and away-side (Ac/) ^ tt) [9[. The combined mea- 
surements of the B decay contribution and of the nu- 
clear modification factor (Raa) for heavy- flavor decay 
electrons in Au-|-Au collisions constrain the value of the 
Raa for electrons from B meson decays. 

The p -\~ p data used in this analysis were taken by the 
STAR experiment [13] during the 2005 and 2006 RHIC 
runs. The main detectors for this analysis are the Time 
Projection Chamber (TPC) and the Barrel Electromag- 
netic Calorimeter (BEMC). The BEMC has a Shower 
Maximum Detector (SMD): proportional gas chambers 
with strip readout at a depth of ~ 5 radiation lengths 
{Xo) designed to measure shower shapes and positions. 
The acceptance for electrons in pseudorapidity and az- 



imuth is \r]\ < 0.7 (0 < 77 < 0.7 in the 2005 run) and 
< (/) < 2tt. The BEMC also serves as a trigger detec- 
tor for high Pt electrons or photons, where single-tower 
transverse energy thresholds of 2.6 GeV and 5.4 GeV 
were used. The total sampled luminosity was 11.3 pb~^ 
(0.65 pb-i) for the 5.4 GeV (2.6 GeV) trigger threshold. 
We used triggered events with primary vertices located 
within 35 cm of the TPC's geometrical center along the 
beam direction. 

Electrons were identified by measuring ionization en- 
ergy loss (dE/dx) and track momentum (p) from TPC, 
the energy (E) deposition in the BEMC, and the shower 
profile in the SMD. A significant fraction of the hadron 
background was rejected by selecting tracks with a mea- 
sured dE/dx in the TPC between —1 and -1-3 standard 
deviations from the expected mean dE/dx for electrons. 
Based on calibrations of the SMD response to electrons 
and hadrons, tracks whose shower projection occupies 
more than 1 strip in both (j) ^-nd 77 SMD planes were 
selected as electron candidates. We required the energy- 
to-momentum ratio to be in the range 0.3 < p/ E < 1.5. 
The hadron contamination in the electron sample after 
applying these cuts is ~ 2% up to 5 GeV/c, increasing 
to - 10% at 9 GeV/c. 

The electron sample has two components: (1) non- 
photonic electrons and (2) photonic electrons - those from 
photon conversion in the detector material between the 
interaction point and the TPC and Dalitz decays, mainly 
from tt". Photonic electrons were identified by pairing 
electrons with oppositely charged partner tracks, deter- 
mining the conversion or decay vertex, and calculating 
the invariant mass of the e^e~ pair Mg+g- jll| . To 
improve the invariant mass resolution, the so-called 2- 
D invariant mass was calculated using only pT and pz, 
which is equivalent to setting the opening angle in the 
transverse plane to zero [ll[. Monte Carlo simulations 
indicate that the cut of 0.1 GeV/c^ removes almost all 
photon conversion candidates for which the decay part- 
ner is reconstructed in the TPC. The efficiencies for pho- 
tonic electron reconstruction (eg ) range from 65% at 3.0 
GeV/c to 80% at 8.0 GeV/c, as determined from GEANT 
simulations. For the enon7-^ analysis, we first removed 
the electrons that have an opposite-sign partner such that 
Mg+g- < 0.1 GeV/c^ from the inclusive electron sample. 
The remaining electrons form the "semi-inclusive" elec- 
tron sample. The non-photonic electron yields can be 
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FIG. 1: (color online) Distributions of the azimuthal angle 
between non-photonic electrons and charged hadrons normal- 
ized per non-photonic electron trigger. The trigger electron 
has (top) 2.5 < pt < 3.5 GeV/c and (bottom) 5.5 < pr < 
6.5 GeV/c. The curves represent PYTHIA calculations for D 
(dotted curve) and B (dashed curve) decays. The fit result is 
shown as the black solid curve. 
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Np . is the number of semi-inclusive electrons. TVpnotrcco 
represents the number of photonic electrons which are not 
reconstructed by the invariant mass method and is de- 
fined as: (l/ee^-l)(A^e„„nk„--^enke)- A^enke IS the uumbcr 
of non-photonic electrons that were rejected by the con- 
version cuts because they happened to form a pair with a 
random track which is determined using like-sign pairs. 
Nfi is the remaining background from hadron contami- 
nation in the electron sample. Other weak decay contri- 
butions such as Ke3 are negligible due to their long ct, 
and charmed baryons (mostly Ac) is expected to be very 
small contribution since the baryon yield is small com- 
pared to the meson yield (Ac/D° ~ 0.1 in PYTHIA) and 
the branching ratio for semi-leptonic decays is smaller 
for baryons than mesons. The Enon-y-h azimuthal distri- 
butions were calculated as 
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where each term is normalised to be per non-photonic 
electron trigger. Each angle-difference distribution on 
the right-hand side of Eq. ^ was experimentally de- 
termined. The distribution dN'^i '''' /d{A(j)) was con- 
structed from dN'^1 ~^ /d{A4>) by removing the conver- 
sion partner to account for the fact that the partner elec- 
tron is not reconstructed. 

Figure [T] shows dA^'^"°""'"''/d(A0) per trigger for non- 
photonic electrons for two different trigger pT selections. 



Associated particles were required to have pr > 0.3 
GcV/c and \r]\ < 1.05. The dotted (dashed) line in the 
figure represents a PYTHIA version 6.22 calculation of 
the azimuthal correlations between electrons from D (B) 
meson decay and charged hadrons {feo{A4>), fegiAcj))) 
12{ . PYTHIA was tuned to reproduce the shapes of pr 
distributions for D mesons measured by STAR [1^, [l3| . 
The PYTHIA calculation shows that the near-side peak 
for fes{A(j)) is broader than that for fe^{A(j)). These 
shapes are dominated by decay kinematics. The frag- 
mentation function does not affect the shape in a signifi- 
cant way. The fraction of non-photonic electrons from B 
meson decay can be determined by fitting the near-side 
distribution function (|A0| < 1.5): 
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= refes (A0) + (1 - re) fen (A0), (3) 



where vb is the ratio of electrons from B meson decay to 
the total non- photonic electron yield, vb = Neg/{Neg + 
Ne,) = NeJNe_.,. 

An independent measurement of vb was performed us- 
ing Cnon^y-D^ Correlations. D^ mesons were reconstructed 
via their hadronic decay D° — !• K~Tr~^ {B = 3.89%) by 
calculating the invariant mass of all oppositely charged 
TPC tracks in the same event. In this analysis, only 
events with a non-photonic electron trigger were used 
for D^ reconstruction. Furthermore, the kaon candidates 
were required to have the same charge sign as the non- 
photonic electrons [9|. The combinatorial background 
of random pairs was evaluated by combining all charged 
tracks with the same charge sign from the same event. 
The requirement of a non-photonic electron trigger sup- 
presses the combinatorial background, yielding a signal 
(S')-to-background {B) ratio of about 14% and a signal 
significance {S/yJS + B) of - 4.6. 

Figure [2] (a) shows the background subtracted ir-K in- 
variant mass distribution. The peak position and width 
were determined using a Gaussian fit to the data. The 
Kit invariant mass distribution was obtained for differ- 
ent A0 bins with respect to the trigger electron, and the 
yield of the associated D^ mesons was taken as the area 
underneath the Gaussian fit to the signal. Figure [5] (b) 
shows the azimuthal correlation of e^o^^-D'^ , which ex- 
hibits near- and away-side correlation peaks with similar 
yields. The results are fitted with the correlation func- 
tions for charm and bottom production from PYTHIA 
and MCQNLO simulations having the relative B contri- 
bution as a free parameter [9|. The observed away-side 
correlation peak can be attributed to prompt charm pair 
production (^ 75%) and B decays (^ 25%), whereas 
the contributions to the near-side peak are mainly from 
B decays. We determined tb by fitting the measured 
e„on7-£'° correlation with PYTHIA and MCQNLO and 
used the average of the two fits for the final value. 

Figure 13] shows r^ — Neg/{Neg +Nen) extracted from 
enon7-^ correlations (filled circles) and Cnon-y-D'^ corrc- 
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FIG. 2: (a) Background-subtracted invariant mass distribu- 
tion of Ktt pairs requiring at least one non-photonic electron 
trigger in the event. The solid line is a Gaussian fit to the data 
near the peak region, (b) Distribution of the azimuthal an- 
gle between non-photonic electron (positron) trigger particles 
and D^{D ). The solid (dashed) line is a fit of the correlation 
function from PYTHIA (MCQNLO) simulations to the data 
points. 



lations (open circle) as a function of px- The verti- 
cal lines represent the statistical errors and the system- 
atic uncertainties are shown as brackets. The system- 
atic uncertainties due to electron identification (~ 7% 
), photonic electron rejection (^ 6%), the fit range 
(~ 10%) and the normalization of the azimuthal distri- 
bution (- 10%), PYTHIA and MCQNLO predictions for 
e-non-y ~ D^ (^ 5%), and the D^ signal extraction were 
estimated by varying the associated cut parameters and 
adding the individual contributions in quadrature, tb in- 
creases with electron pr and reaches approximately 0.5 
[Neg/Ne^ ^ 1) around pt = 5 GeV/c. vb from the 
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D° correlation measurement at px ^ 5.5 GeV/c is 
consistent with tb from 6,10117-^ correlations. The curve 
in the figure is tb from a FONLL pQCD calculation in- 
cluding theoretical uncertainties |14l |. Similar ratios at 
2 < px < 7 GeV/c using a different method have also 
been reported [15[. J/^p di-electron decays can also con- 
tribute to non-photonic electrons and STAR measure- 
ment of J/ip at high pt indicates that J/ip decays could 
contribute nearly 10 % around pr 5 GeV/c. The esti- 
mated effect of the electrons from J/ip decays on tb is 
a few percent, much smaller than the current statisti- 
cal and systematic uncertainties, and no correction was 
applied to our data. 

Next, we explore the implications of the measured r^ 
for the Raa of electrons from B meson decay in heavy 
ion collisions. The Raa for heavy flavor non-photonic 
electrons {R^^ ) is given by 



RTA^a-rB)RTA + rBR'A% 



(4) 



where i?2^ {^aa) i^ ^^e Raa for electrons from D (B) 
mesons. From Eq. (|3]), i?^^ ^^'^ -^AA ^^^ related by 
the B decay contribution to the non-photonic electron 
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FIG. 3: Transverse momentum dependence of the relative 
contribution from B mesons (rs) to the non-photonic electron 
yields. Error bars are statistical and brackets are systematic 
uncertainties. The solid curve is the FONLL calculation [l4l ]. 
Theoretical uncertainties are indicated by the dashed curves. 



yields (r^) in p + p collisions. We have taken the i?^^ 
measurement from PHENIX [H [l^ and fit the R% 
above px > 5 GeV/c to a constant value and obtained: 
Raa = 0.167t°:°ti (stat)+°;°^;^ (syst) ±0.0117 (norm), 
where the statistical and systematic errors are evalu- 
ated from weighted average over these pr > 5 GeV/c 
points. We also calculate the weighted mean rs value 
for Pt > 5 GeV/c including statistical and systematic er- 
rors from our measurement: r^ = 0.54 ± 0.0349 (sta.) ± 
0.0666 (sys.). Then using Eq. 4 we calculate a likelihood 
distribution for -R^^ as a function of i?^^ and the re- 
sults are shown in Fig. 4. The most probable values for 
the i?^'^ and i?^^ correlation are shown by the line with 
open circles and the 90% Confidence Limit curves are 
represented by dashed lines. This result indicates that 
B meson yields are suppressed at high pr in heavy ion 
collisions presumably due to energy loss of the b quark 
in the dense medium [5| or of the heavy flavor hadrons 
due to dissociation [6| or elastic scattering [7| . Our con- 
clusion does not change if we use the Raa measurement 
from [l| and ignore the J/tjj feeddown contributions. 



For comparison, we also show model calculations in 
Fig. 01 Model I includes radiative energy loss via a 
few hard scatterings with initial gluon density dNg/dy 
— 1000 y. Model n includes cold nuclear matter ef- 
fects, partonic energy loss and collisional dissociation [6|. 
Model in assumes a large elastic scattering cross section 
associated with resonance states of D and B mesons in 
the QGP [71 . The model contours in Fig. |4]are calculated 
from the pr dependences of Raa for D and B decay in 
the interval 5 < pt ^ 9 GeV/c. For model I and II, the 
uncertainties are also taken into account. The experi- 
mental results are consistent with models II and IH but 
are incompatible with model I. Recently AdS/CFT the- 



< 1 

0) < 

CCo.g 



0.8 
0.7 
0.6 

0.5 
0.4 
0.3 
0.2 
0.1 
0, 




I I l:Djordjevic 
I I IliSharma 
— Ilhvan Hees 

*R™™'<(0-10%) 
CL 90 % 



R 



AA 



cncc Grid consortium for providing resources and sup- 
port. This work was supported in part by the Offices 
of NP and HEP within the U.S. DOE Office of Science, 
the U.S. NSF, the Sloan Foundation, the DFG cluster of 
excellence 'Origin and Structure of the Universe' of Ger- 
many, GNRS/IN2P3, STFG and EPSRC of the United 
Kingdom, FAPESP CNPq of Brazil, Ministry of Ed. and 
Sci. of the Russian Federation, NNSFC, CAS, MoST, 
and MoE of China, GA and MSMT of the Czech Repub- 
lic, FOM and NWO of the Netherlands, DAE, DST, and 
CSIR of India, Polish Ministry of Sci. and Higher Ed., 
Korea Research Foundation, Ministry of Sci., Ed. and 
Sports of the Rep. Of Croatia, Russian Ministry of Sci. 
and Tech, and RosAtom of Russia. 



FIG. 4: (color online) Confidence level contours for nuclear 
modification factor Raa for electrons from D (R'^a) and B 
( -R^^) meson decays and determined by combining the Raa 
results and the rs measurement for pr > 5 GeV/c. Three 
different models of Raa for D and B are described in the 
text. 



ory has also been used to calculate the heavy quark en- 
ergy loss in a str ongly c oupled quark-gluon plasma mat- 
ter, for example 



18l Il9j . The theory also predicts strong 



suppressions for charm and bottom, and the ratio of the 
nuclear modification factors is proposed to differentiate 
AdS/CFT calculation from others j20| . 

In summary, we measured the relative contribution 
from B decays to the non-photonic electron production in 
p + p collisions at -^/i = 200 GeV by using azimuthal cor- 
relations between non-photonic electrons and hadrons {h, 
D°). Our result indicates that the B decay contribution 
increases with px and is comparable to the contribution 
from D meson decay atpT>5 GeV/c. Our measurement 
is consistent with the FONLL calculation. The ratio of 
N(.g/{Ne^ +Neg ) combined with the large suppression of 
non-photonic electrons indicates that Raa for electrons 
from B hadron decays is significantly smaller than unity 
and therefore B meson production is suppressed at high 
Pt in heavy ion collisions. The constraint on Raa ^^'^ 
R'Aji will help to differentiate theoretical model calcula- 
tions for heavy quark energy loss in the dense medium. 
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